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In the usual parameter regime of accelerator physics, particle ensembles can be treated as classical. 
If we approach a regime where CxeyCe ~ ■^particiesAcomptom however, the granular structure of 
quantum-mechanical phase space becomes a concern. In particular, we have to consider the Pauli 
exclusion principle, which will limit the minimum achievable emittance for a beam of fermions. We 
calculate these lowest emittances for the cases of bunched and coasting beams at zero temperature 
and their first-order change rate at finite temperature. The self-field of the particle beam will, 
analogous to a space-charge tune depression, lead to a decrease of the quantum-mechanical state 
density. We calculate the tunes and the increased emittance for this case. 



I. INTRODUCTION 

When a relativistic particle travels in a focusing dominated beamline, it can be viewed, in its rest frame, as a 
particle in a harmonic oscillator potential. Emission of dipole radiation will lead to transverse energy loss; quantum- 
mechanically, the particle will drop down one energy level by emission of a photon. This mechanism will, in the 
absence of competing heating mechanisms, ultimately lead to a cooling down of the particle's transverse action to its 
quantum- mechanical limit of Jmin — ^- This result has been established in 0|, however, only single-particle dynamics 
is considered there. 

In the multi-particle case, some care has to be taken if the particles in question are fermions, as Pauli's exclusion 
principle has to be considered when constructing the ground state; it has been claimed elsewhere that this will limit 
the admissible minimum emittance to values not too far out of the reach of today's technologies. 

In this paper, the emittances for a relativistic fermion bunch and a coasting beam are calculated for temperatures 
close to zero. We will show that the quantum-mechanical limit for the transverse emittances for typical configurations 
are vastly smaller than values reached today. 

II. CALCULATING COLLECTIVE QUANTITIES: THE WIGNER TRANSFORMATION 

In this paper, we want to calculate the emittance of a multi-particle system. The classical emittance is defined in 
terms of classical expectation values, i.e. particle averages, as = {q'^)^ " (w)ci- ^ quantum version of this 

and similar quantities can be obtained by interpreting the Wigner transform of the state of the quantum system in 
question as a classical probability distribution. 

Consider an A^-particle system in a pure state l^*) = P \ni) ■ ■ ■ |n„), where P is the (anti)symmetrizing operator, 
the particles in question being (fermions) bosons. Now, we can create a pseudo-classical probability distribution 
p{QijP1i • • • ; QNiPn) by using the Wigner transform of the product wavefunction. It is easy to see that all expectation 
values of this phase-space distribution are the same as the ones of the Wigner transform of a one-particle statistical 
matrix with equal weights in all states. We can thus conclude that the appropriate quantum version of the classical 



X. 

H , expectation values above is obtained by replacing the classical average {f{p,q)) by ( (/(Pop, 9op))q,„ ) , i-C, the 
quantum-mechanical expectation value for the corresponding operator function, averaged over the occupied states. 



III. DYNAMICS 



Consider an ultra-relativistic particle beam in a circular accelerator. Neglecting higher-order effects, the Hamilto- 
nian can be written as a quadratic form in the usual phase-space coordinates x, x' , y, y' , a, S. 
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However, this Hamiltonian is not appropriate for quantization, as energy and time have switched roles. Thus, we 
use the Hamiltonian of the system in the beam's frame of reference, which can be obtained by a series of canonical 
transformation from the lab frame 

where 7 is the relativistic factor, Kx.y the (possibly local) focusing strengths (in the case of magnetic quadrupoles, one 
has Kx = ~i^y) E^nd the external electric potential. The directions x,y,z are radial, transversal, and tangential, 
respectively. Note that we use units with h — c — ks — mo throughout, so all quantities are expressed in powers of 
the Compton length of the particles in question. 

The longitudinal part of the Hamiltonian depends on the physical setup. The particles might either be confined 
by the nearly harmonic potential of the RF bucket, or we have case of a coasting beam, where the only constraints 
imposed on the longitudinal motion are the ones due to the periodicity of the problem. In the sequel, we will consider 
both cases. 



IV. ANISOTROPIC OSCILLATOR 



Let us assume that the longitudinal motion is determined by the presence of an RF bucket. We can approximate 
the potential $ by expanding it to 2nd order in z. For reasons of simplicity, we only take into account the O(z^) 
term, i. e., we assume that the particle is on the orbit and is not losing energy. 

For a bunched beam with dimensions cr^, as, we have, in the orbit frame, toijai — as. 

The longitudinal and radial parts of the hamiltonian (^ have the form 

2 2 22 

H = 1^Y ~^ — '"^^^^ ■ ^ ' 

i=l 

Its eigenfrequencies are determined by the equation 

+ X^{loI+loI)+loI{^? -LoD^Q , (3) 



which leads to stable motion for 



{ujI + ljI) > 4w|(^^ ~ujI) >0 



The first condition will always be fulfilled for realistic machines. The second one corresponds to the machine being 
below or above transition: if the second factor changes sign, the eigenfrequency can be mad real again by flipping the 
sign of ^'rf- However, the absolute sign of both the kinetic and the potential term will change, leading to the 
(for purposes of constructing the quantum-mechanical ground state) pathological case of a hamiltonian not limited 
from below. In the sequel, we will assume the machine is below transition. 

Thus, expanding (|^) to first non-trivial order in the canonical coordinates and applying the canonical transformation 
removing the mixed term in (^, we obtain the Hamiltonian of a 3-dimensional harmonic oscillator with corrected 
frequencies given by (|3|); the ground state is characterized by the occupation numbers njj e {0, 1} where '^d = ^f 
the oscillator levels. For sake of generality, we consider the case of d dimensions. The ground state for a given particle 
number can be constructed by successively filling states with the lowest energy (we disregard spin here, which can be 
easily reintroduced by replacing N — > 2N in the final formulae) . 

In ^-space, the Fermi sea is just a unit d-simplex, in rii-space, a d-simplex with axes of length . . . Thus, 
the particle number for a ground state filled up to the Fermi energy cp, where we have disregardecf the zero- mode 
energy i uji of the oscillator, 

d 

nd\ ' 

the volume of an unit d-simplex being ^ and ft = uji ■ ■ ■ uJd- 

The energy in the ith degree of freedom in that case is given by a sum over the d-simplex. Replacing all sums by 
integrals, we have 



= IT / ■ ■ ■ / EFgidqi ■■■dqd = 



n Jo Ja Jo ^ ^ " n{d + iy. d + i 
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We calculate the emittance by using its statistical definition = (p^) <f — (pq)^ and replacing the classical averages 
by the ones discussed in II. 

For a harmonic oscillator, we have (p^)qjjj ('Z^)q,„ — ^^d (m)qi„ — 0, so: 

-/ \ - 

and the total phase-space volume 



Thus, the projected emittances scale as TVi, as one would naively assume. Furthermore, due to the occurrence of 
the geometric mean of the frequencies in (^), the projected emittance in one dimension can be lowered by shallowing 
the potential in the other dimensions. 

Note that a similar approach has been chosen elsewhere; [|^] gives an estimate for £min from a similar reasoning, 
but ends up (due to a miscounting of the states) with a scaling different from our result. 



'Nd\ 



(4) 



n 



= £, = 



Ndl 



V. MIXED CASE: LONGITUDINALLY FREE PARTICLES 



So far, we have assumed an anisotropic oscillator. But given the case of a particle moving freely longitudinally, the 
energy content of that degree of freedom will be given by the square of the (angular) momentum. (We might consider 
the boundary conditions imposed by a periodic box instead of a circular arrangement.) We treat the general case, 
i. e., a hamiltonian 



LUifli 



The emittance in each of d new degrees of freedom is given by {p^)q^ — = 2£'„ and {<i'^)q^ — where D 
is the length of the enclosing potential, and (M)Qm ~ 0' 



D 



\ 



\ V|„) 



= JE: 



6 V \ '/|">\/12<D- 

Rescaling the integration range to a unit sphere and unit simplex, we get 



N = 



dqdq 



i70 J (i-sphcrc J rf-simplcx 

We can readily write down the averaged values of the energy in the different degrees of freedom: 



(5) 



£^7 



and 



{Ei. 



= if) = 



-+Mq 



2d + d + 2 



+d-l 



dq 



J^{l-q)i+'i~hlq 2d + d + 2 



(6) 



uJi d + 2d+2 



(7) 
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and 



and the product emittance is 



12(d + 2d + 2) 



dr,d '^+2 / TT \ 4 ^ [9 
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Of course, by putting d = 0, we regain the formulae for the bunched-beam case. 

Putting in a real coasting-beam ring, we can express the orbit-frame frequencies by the tunejj]: uix ~ — 
where L is the length of the ring. The longitudinal momentum is quantized in units of , so = ^^j^ and 



uJx 7 1/ 256'y Li^ 

This means the transverse degrees of freedom begin to exhibit a non- minimal emittance for N > Nq, where 



iVo ~ 78.0. 



1/7 L 



A, 



Compton 



in usual units, and the transverse emittances will grow cx N's for higher particle numbers. 

For parameters one could consider for a real focusing-dominated, below-transition ring (7 = 10, v — 100, L — 27rm), 
we get A'o « 1.0 • 10^" particles; i.e., the beam will have, for realistic particle numbers, an emittance close to its 
minimum value. (Acompton in customary units.) 



VI. SELF-CONSISTENT TUNE SHIFT FOR THE FERMI CONDENSATE 



In our construction, we tacitly assume that the particle-particle interaction does not modify the particle content of 
the ground state. This corresponds precisely to the notion of a Fermi liquid (in our case, a highly anisotropic one), 
in which the free particle spectrum smoothly deforms into the quasi-particle spectrum of equal particle content when 
the interaction is switched on adiabatically. This naive assumption of the existence of a Fermi surface may break 
down if we take into account particle-particle interactions. 

For the case of the system being below transition, we can make the following semi-quantitative argument for the 
existence of a Fermi liquid: In a "mean-field calculation" , we estimate the effective transverse focusing strength Woff 
to be the sum of the external focusing and a space-charge tune depression due to a circular beam of radius -y/ (x^) : 

2 2 2 2 JVC 

^eff = '^ext - '^SC = ^oxt " „ / x 2 ' 

27L [x) 

In the foUwing, we consider a coasting beam. We have (^x'^') cu'^g = {Ex), where {Ex) is given by (||). We use (||) 
with uj — LUcft to eliminate ep and obtain a consistency condition for ujcS- 



where we have introduced 



'J-'ext 



^ / Whrt / ^cff \ 



-4/5' 



V '^ext / \ l^ext / 

{/2375r,V7£)jV3 

IS-TT^Acompton 



(8) 



(9) 



(||) has a solution Wpff < i^oxt for for all N. This corrected frequency has to be substituted in the equations in the 
previous section, leading to corrected values for the emittances. 
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FIG. 1: Effective focusing strength of a coasting beam in its ground state; ujint ~ e^N /{■yL) 
For typical setups, Wint ^ Wcxt 

(using the ring parameters from above and N — 10^*^ particles, we end up with cjint/woxt ~ 2.14-10^ and Ci, corrected ~ 
18.5£i), so ~ (w^) ■ ^^'^ scaling of the transverse emittances is oc ui^^\ with ujcs « N^i , so we end up 

with a corrected scaling for the case of large particle numbers: oc instead of cx Ni from (0). 

This modified scaling law means that the transverse phase space volume scales as SxSy — 0{N): adding a new 
particle to the coasting beam's ground state will lower the transverse frequencies, thus increasing the state density, so 
a new site becomes available in transverse phase space, while the longitudinal phase space retains its original volume. 
In particluar, this mes-iis ths-t kinetic cind potentiRl energy cilws-ys will b^la^nce out, cp ~ ^kin ~ ^poti 

keeping the 

system at the boundary of a "weakly non-ideal" and a "strongly non-ideal" gas. Thus, the simple criterion efei„ <C epot 
for quantum Wigner crystallization is not applicable for any external focusing strength. This property is peculiar to 
the external oscillator potential and is not present for other confining potentials. 

VII. FINITE TEMPERATURE 

The above considerations were for the case of zero temperature. To generalize to finite temperatures, we follow 
the usual prescription and introduce a chemical potential. The quantity we want to calculate is the logarithm of the 
partition function of the grand-canonical ensemble: 

oo 

logZ = log ^ e-'^S*"^(S'=i"'=(''= + 5)-^) = ^ log(l + e-''(^'"^(*'' + 5)-f)) 

Again, we transform the sum into an integral. With (|^), we integrate over energy 

log Z = i ^ r log f 1 e-'3(^-^)) E^+i-^dE . 



For small temperatures, integrals of this style can be done by integrating by parts and then be approximated using a 
Sommerfeld expansion. We find 



\ogZ 



d+l 



As 



_l d_ 



\ogZ 



d 



^logZ 



we can write the temperature-dependent contributions to the emittances 



+ 1 



6/32 



(10) 



where we have used the zero-temperature ep as chemical potential, which is correct in this order of /?. Thus, the 
emittances grow quadratically in the temperature, with the scale being set by e^. For the parameters from above and 
N = 10^" particles, one finds, from (|lO|), a transverse emittance-doubling temperature of T = 0.03K. 
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